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We investigate numerically with Hylleraas coordinates the frequency dependence of the carrier- 
envelope phase (CEP) effect on bound-bound transitions of helium induced by an ultrashort laser 
pulse of few cycles. We find that the CEP effect is very sensitive to the carrier frequency of the laser 
pulse, occurring regularly even at far-off resonance frequencies. By analyzing a two-level model, we 
find that the CEP effect can be attributed to the quantum interference between neighboring multi- 
photon transition pathways, which is made possible by the broadened spectrum of the ultrashort 
laser pulse. A general picture is developed along this line to understand the sensitivity of the CEP 
effect to laser's carrier frequency. Multi-level influence on the CEP effect is also discussed. 

PACS numbers: 32.80.Rm, 42.50.Hz, 32.80.Qk 



INTRODUCTION 

For an ultrashort laser pulse that lasts for only a 
few cycles, its carrier-envelope phase (CEP) can dra- 
matically affect the yield of matter-laser interactionhj, 
leading to CEP dependence of electron ionization [2j-|6[ 
and harmonic-photon emission 

0-Eoj]. Recently, with 
the rapid development of the laser technology, CEP has 
become a new way to control the dynamic process of 
matter-laser interaction [9]. It has been demonstrated 
that the CEP effect for an intense laser pulse can be 
measured by comparing the photoelectron yields in two 
opposite directions along the laser's electric field 0. 

More recently, the CEP effect on bound-bound transi- 
tion of an atom has been investigated theoretically fill — 
13 1 and observed experimentally [14]. Roudnev and Esry 
have presented a general framework for understanding 
the CEP effect using the Floquet theory Li et al. 

have demonstrated that an experimentally observed CEP 
effect can be attributed to the interference between one- 
and three-photon transition pathways LJ]. The study 



by Nakajima and Watanabe suggests that the CEP ef- 
fect can occur as the laser's carrier frequency is far off- 
resonance 



Hi- 



In this paper, we use Hylleraas coordinates to study 
numerically how the CEP effect on bound-bound transi- 
tions of helium changes as a function of the carrier fre- 
quency of an ultrashort laser. Our numerical results show 
that the CEP effect depends sensitively on the carrier fre- 
quency even when it is far off-resonance. The essential 
physics implied in these numerical results can be well re- 
vealed by a two-level model. For a two-level system, when 
the pulse duration is long, quantum transitions peak 



at well-separated multi-photon resonant frequencies. As 
the pulse duration decreases to less than three laser cy- 
cles, for example, the widths of such transition peaks get 
significantly broadened, and eventually two broadened 
neighboring transition peaks can cross with each other. 
As a result, two different multi-photon transitions can 
both contribute significantly to the total transition am- 
plitude and hence interfere with each other. We show 
that a large CEP effect occurs exactly at these cross- 
ings and can hence be understood as a quantum interfer- 
ence effect [17j . A general and simple picture developed 
along this line enables us to clearly answer the following 
questions: (1) Why is the CEP effect on bound-bound 
transitions sensitive to a carrier frequency that is far-off 
resonance? (2) Which multi-photon transition pathways 
can interfere and lead to the CEP effect? 



NUMERICAL METHOD TOWARDS THE CEP 
EFFECT IN HELIUM 

Atomic units are used throughout unless specified oth- 
erwise. To study a photoexcitation process of a ground 
state helium atom in a linearly polarized ultrashort laser 
pulse, we solve a time-dependent Schrodinger equation 
= [Ho + Hi(*)]|#(t)). The field-free part of the 
Hamiltonian reads 
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where ri and r 2 are the coordinates of the two electrons 
measured from the nucleus located at the origin. The 
light-atom interaction part of the Hamiltonian is 

Hi = (n + r 2 ) • e(t) = -(ri + r 2 ) • dA(t)/8t , (2) 



where the vector potential of the field is given by 0, EH 

A(t) = sA exp(-a 2 t 2 ) sin(wt + 4>) /w, (3) 

with e being the linear polarization vector, u the carrier 
frequency, and <fi the CEP parameter. The full width 
at half maximum (FWHM) of the pulse duration is r = 
2v / m~2/a. 

The wave function of helium is expanded in 

terms of the field- free eigenvectors \ip n ): = 
2 n o„(t)e" ,B "'|^) , where H \tp n ) = E n \tp n ). Then we 
numerically solve the following set of equations for the 
amplitudes a n (t) 
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where S mn = # m - £? n , and H^ m (t) = (ip n \H.i\ip m ) . 
Initially, the helium is in the ground state. The basis set 
is constructed in Hylleraas coordinates [l6[ 



H'n) = Cnm\<l>n 
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Figure 1: (color online) The CEP parameter M. [see Eq.©] 
versus laser's carrier frequency for the 2 1 P state (a) and the 
3 D state (b) of helium. Solid squares and solid circles are 
numerical results, (green) dashed line is for a two-level model, 
(blue) dot-dashed lines are for a three-level model, and (red) 
solid lines are for a four-level model. 



where 



\<Pn 



In the above, yPf 1 (ri, r 2 ) is the vector coupled product 
of spherical harmonics for the two electrons. In order 
to obtain convergent results, 874 Hylleraas functions are 
used with the total angular momentum quantum number 
L ranging from to 4. 

Our numerical results show that the population of an 
excited state of helium can strongly depend on the CEP 
parameter <fr at certain carrier frequencies. To quantify 
this CEP effect, we introduce the parameter 11 1 



M = 



-P( 







(-P(0max) + P(&ni„))/2 ' 



(5) 



where P(</> m ax) and P(^> m in) are, respectively, the maxi- 
mum and minimum populations for a given excited state. 
A large value of M. means a strong CEP effect. We nu- 
merically calculate M. for the 2 1 P and 3 1 !) states. In 
our computation, the FWHM pulse duration is one laser 
cycle and the peak intensity of the laser is 10 13 W/cm 2 . 

Figure Q] presents how A4 changes with the laser car- 
rier frequency u. Results for the 2 1 P state (squares) 
are shown in panel (a) and those for the 3 1 !) state (cir- 
cles) are shown in panel (b). It is clear from the figure 
that there are many carrier frequency windows in which 
M. peaks at about its maximum value of 2, indicating 
a strong CEP effect on the transition probabilities from 
the ground state to the two excited states. These nu- 
merical details qualitatively agree with those reported in 
Ref. 14j, where the authors changed the energy differ- 



One important feature in Fig. Q] is that the CEP ef- 
fect is very sensitive to the carrier frequency even it is 
far off-resonance as compared with the energy difference 
between the two bound states. This implies that far off- 
resonance is not a sufficient condition for a large CEP 
effect. In the following, we attempt to answer the follow- 
ing question: what is the physics underlying these narrow 
frequency windows where the parameter Ai peaks? 

We have compared our full numerical results in 
Fig. [1] (a) to the results of a two-level model (to be elab- 
orated below) . As clearly shown in the figure, except the 
peak positions are shifted towards lower frequencies, the 
two-level results (the (green) dashed curve in Fig. D](a)) 
can embody the main features in the full numerical re- 
sults for u > 0.1. This indicates that a simple two- level 
model is sufficient to reflect the essential physics behind 
these peaks. In the following, we first study a two-level 
model, and then investigate multi-level influence on the 
CEP effect. 



TWO-LEVEL MODEL AND INTERFERENCE 
BETWEEN NEIGHBORING TRANSITION 
PATHWAYS 

Consider a two-level system in a pulsed laser field. The 
ground and the first excited states of helium are chosen 
as the two levels, denoted as |0) and |1). According to 
Eq. (J4j) , the amplitudes of these two states obey the fol- 
lowing equations 



ence between two bound states instead. 
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where \ijk = Hkj represents the transition dipole moment 
between two quantum states \j) and |fc) and fjk(t) = 
e(t) exp(iAEj k t) with e(t) being laser's electric field and 
AEj k the energy difference between |j) and \k). 

If the system is initially in the ground state, a formal 
solution of a!, after the laser pulse has passed, can be 
written as 

/oo 
dtfi (t) 
-OO 

/CO pt ptl 
/ / f w(t) f 01(h) fwit^dtdhdtz 
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+ ■■■ 

OO 

= Y^T^ 2n+1 \ (7) 

n=0 

where T( 2ll+1 ) represents the (2n + l)-photon quantum 
transition amplitude. For example, is for one-photon 
transition amplitude and is for three-photon transi- 
tion amplitude. 
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j>(2n+i) can jj e decomposed into a sum of (n + 1) terms, 
i.e., 

n 

T (2n+1) K J2 T SlV\ (8) 
3=0 

where 

T^l^ cx 5 [AE 10 - (2j + l)w] exp[-i(2j + 1)4>] . (9) 

This implies that, physically, each term T^ 1 ^ 1 ^ can be 
associated with a quantum transition pathway, which 
carries a phase — (2j + 1)4> and contributes significantly 
to the overall transition amplitude at resonant frequency 
ujj = AEio/(2j + 1). These quantum transition path- 
ways are depicted schematically in Figure [U where the 
roth column shows all the pathways involving (2m — 1) 
photons and the kth row includes all the pathways that 
have a resonant frequency at tuk-i = AE\q / (2k — 1) with 
ro, < 4 and k < 3. 

With the decomposition in Eq.©, it is clear that the 
overall transition amplitude a\ is significantly different 
from zero only at resonance frequencies. At the same 
time, we notice that all the pathways at a given resonance 
frequency ujj carry the same phase —(2j + l)(j). Therefore, 
the CEP </> does not affect the magnitude of a±. This 
clearly explains why there is no CEP effect for a long- 
pulsed laser. 
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Figure 2: Quantum transition pathways between two bound 
states of an atom. For the pathways in the same column, they 
involve the same number of photons, e.g., the third column is 
for five-photon pathways. For the pathways in the same row, 
they share the same resonance frequency, e.g., the pathways in 
the second row have the resonance frequency A_Bio/3, where 
AE10 is the energy spacing between the two bound states. 

To clearly understand the physics behind each ampli- 
tude T ty2n+1 \ we first consider the limit of long pulse 
duration, i.e., a —¥ 0, which is easier to deal with. In 
this limit, after neglecting the negative frequency com- 
ponents [1J|, the (2?t, + l)-photon transition amplitude 



Figure 3: The mag nitude of T (2n+1) with n = (solid curves), 
n = 1 (dashed curves), and n = 2 (dash-dotted curves) for 
pulse duration (a) r = 20, (b) r = 7, (c) t = 3, and (d) r = 1 
laser cycle(s). The arrows in panel (c) indicate the interfer- 
ence between T 3 (3) and t[ 3) (right) and T 3 (3) and T 5 (5) (left); 
the arrows in panel (d) indicate the interference between 
and T 3 (3) (right) and T 3 (3) and T 5 (5) (left). 

The situation becomes very different as the laser pulse 
becomes shorter. It is reasonable to assume that the 
decomposition in Eq.© still holds even for short laser 
pulse. However, as the pulse becomes shorter, T^+i^ ^ s 
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no longer proportional to a delta function as in Eq.©, 
but becomes a function of the laser's frequency which 
peaks at (Jj = A_E 10 /(2j + 1). Moreover, the width of the 

peak associated with each transition path T 2 ^"^ broad- 
ens as the pulse duration gets shorter. Consequently, 
when the pulse is very short, e.g., lasting only for a few 
cycles, the peaks can become so wide that the peaks for 
pathways T 2 j+i^ with different j's can cross with each 
other at a certain far off-resonance frequency. As a re- 
sult, the overall transition amplitude a\ can be regarded 
as an interference between the two pathways, if it is dom- 
inated by this pair of pathways at the crossing point. At 
the same time, we notice that each pathway T^j^ car- 
ries the phase — (2j + 1)0, indicating that the relative 
phase between this pair of pathways depend on the CEP 
</>. These facts mean that the overall transition amplitude 
a\ depends on (f> at the crossing point and is strongly af- 
fected by the CEP. 

The above analysis is confirmed by our detailed numer- 
ical calculations shown in Fig. [3j where we analyze the 
amplitude R = |T^ 2 ™ +1 ^| as a function of the laser's car- 
rier frequency for different pulse durations while the CEP 
is fixed at </> = 0. As clearly demonstrated in Fig. Ela), 
for a laser pulse of 20 cycles, the peaks at resonance fre- 
quencies wo, Ui, and u>2 are narrow and well separated. 
Particularly, in consistent with the above picture for the 
long pulse limit, one finds that the mag nitude of 
has one peak located at ojq, indicating the contribution 
from ; whereas the magnitude of has two peaks 
located at u>o and u)\ , indicating the respective contribu- 
tions from T x (3) and T 3 (3) . Finally the magnitude of T< 5 ) 
has three peaks located at u>q, uj\, and u>2, due to , 

(5) (5) 

Tg , and Tg , respectively. As the pulse is shortened 
to seven cycles, there is no essential change in the over- 
all features: only each peak becomes wider, as demon- 
strated in Fig. EJb) , which shows that the decomposition 
for short pulses in Eq.© is well justified. 

It is a very different situation when the pulse is short- 
ened to 3 cycles. In this case, the peaks become so broad 
that they begin to overlap and cross into each other. As 
indicated by arrows in Fig.[3](c), there are two crossings. 
One occurs at w — 0.380, where the peak associated with 

(3) (3) 

Tf crosses with the peak associated with T 3 . The 
other crossing happens at ui = 0.185, where the peak 

(3) 

associated with T 3 crosses with the peak associated 
with T 5 (5) . At such crossing points, the two dominant 
pathways have the same magnitude, which is essential 
for effective interference and strong CEP effect. Finally, 
when the pulse duration is decreased to one laser cycle in 
Fig. [3](d), the contributions from the pathways and 

(3) 

T 3 cross with each other at u = 0.277 and the contribu- 
tions from T 3 (3) and T 5 (5) cross at w — 0.127. Compared 
with the two-level results in Fig. [TJ it is seen that the 
CEP effect occurs exactly at the frequencies where the 



crossings between the broadened peaks occur. 
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Figure 4: Multi-photon transition contributions to the excited 
state amplitudes r pi 2n+1 ) with n = 0,1,2, as a function of 
laser's carrier frequency. The CEP parameter is given by 
4> = in panels (a) and (b), and <f> = tt/2 in panels (c) and 
(d). 

To see how the CEP affects the transition ampli- 
tude, we next investigate both the magnitudes and the 
phases of T^ 3 \ and by setting T( 2n+1 ) = 

i? T (2n+i) exp (if3 T (2n+i) ) with the pulse duration fixed at 
one laser cycle. The results are shown in Fig. |4|a) and 
(b) for cj> = and Fig. H^c) and (d) for = tt/2. Ar- 
rows in Fig. |4] (a) indicate the two frequencies w = 0.277 
(solid arrows) and w = 0.127 (dashed arrows) at which 
the CEP effect occurs. At w = 0.277, where the curves 
for R T {i) and R T (si cross with each other, we see from 
Fig. H] (b) and (d) that the phase difference 8 T (i) — 9 T ( 3 ) 
is around tt at (f> — and around zero at cj> = tt/2. This 
means that when cj> = tt/2, the contributions to the total 
transition amplitude from and constructively 
interfere, whereas at <j) — they destructively interfere, 
leading to a significant CEP effect at ui = 0.277. There- 
fore, the CEP effect is the result of the ^-dependent in- 
terference between the one-photon contribution and 
the three-photon contribution . Similarly, the strong 
CEP effect at lu = 0.127 can be attributed to the <fi- 
dependent interference between the three-photon ampli- 
tude T^ 3 ) and the five-photon amplitude T^ 5 \ 

Based on this simple picture, it is apparent that the 
CEP effect often appears at a carrier frequency that is 
off resonance: the contributions to the total transition 
amplitude from two different pathways with different j's 
can be comparable only at a frequency that is not equal 
to the resonance frequency u)j. 

In principle, if the pulse duration is short enough, the 
transition amplitudes, which are contributed by any two 
neighboring pathways associated with different resonance 
frequencies Wj's, can interfere with each other. As such, 
one might expect the curve of A4 vs uj to have an infinite 
number of peaks. However, since the transition ampli- 
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tude associated with a multi-photon pathway decreases 
rapidly with the number of photons involved, the actual 
number of clear CEP peaks will be limited by the highest 
order of multi-photon transitions that can have a signifi- 
cant transition amplitude. For example, in our two-level 
model with the laser intensity adopted, the highest multi- 
photon order is five, thus yielding only two peaks in the 
curve of M. vsu; these two peaks arise from the interfer- 
ence between one- and three- photon pathways, as well as 
between three- and five-photon pathways. If the laser in- 
tensity is further increased, then more CEP peaks can be 
expected. Scully et al. [14J presented an example of the 
CEP effect caused by the interference between the one- 
and three-photon pathways, which is just one of many 
possible peaks in our general picture. 
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Figure 5: The frequency versus the laser intensity (a) and 
the pulse duration (b) at which the CEP effect occurs in the 
two-level model. 

We next describe some computational details regard- 
ing how the laser intensity and the laser pulse duration 
change the CEP effect. Figure EJa) depicts how the two 
CEP frequencies (i.e., the carrier frequencies at which the 
M. can reach a large value) vary with the laser intensity, 
where the pulse duration is fixed at r = 1 laser cycle. 
Figure [5jb) shows the impact of the pulse duration, with 
the laser intensity fixed at 10 13 W/cm 2 . It is seen from 
Fig.E^a) that the CEP frequency increases slightly as the 
intensity increases. This is because the magnitude of the 
multi-photon transition contribution yO^ 1 ) is propor- 
tional to /( 2n + 1 )/ 2 , with I being the laser intensity. As 
such, the amplitudes of higher order multi-photon tran- 
sitions increase with the laser intensity much faster than 
those of lower order ones, resulting in that the crossing 
of neighboring peaks associated with different transition 
pathways occurs at higher frequencies as the laser inten- 
sity increases. Turning to the pulse duration dependence 
shown in Fig. Efb), it is observed that the CEP frequen- 
cies are also slightly blue shifted as the pulse duration in- 
creases, indicating that the crossing between neighboring 
peaks occurs at slightly higher frequencies as the pulse 



duration increases. When the pulse duration is larger 
than about four laser cycles, the CEP effect disappears 
for I = 10 13 W/cm 2 , implying the vanishing of cross- 
ing between any two neighboring multi-photon transition 
pathways. 
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Figure 6: The magnitude R of , as well as the magnitudes 
of the transition amplitudes formed respectively by the sub- 
path 0-1-0-1 (dashed line) and by the sub-path 0-1-2-1 (solid 
line) associated with T' 3 ', as a function of the laser's carrier 
frequency. The results are obtained using a three-level model. 
In panel (a) the CEP parameter is given by <j> — and in panel 
(b) the CEP parameter is given by <f> — n/2. 



MULTI-LEVEL EFFECT 

For a real atom, laser-atom interaction induces transi- 
tions among many states. As such, we need to consider 
the influence of multi-level transitions on the final pop- 
ulation of the state of interest. As shown in Fig. Q] (a), 
the results of the two-level model disagree quantitatively 
with the full numerical calculations: the two peaks at 
u! = 0.277 and uj — 0.127 appeared in the two-level model 
are red-shifted as compared to the peaks at uj = 0.315 
and uj = 0.15 found in the full numerical approach. One 
is then motivated to examine a three-level model that 
incorporates |0), |1), and \2) = 3 1 !? states of helium. 
Interestingly, the first peak found in this 3-level model 
agrees well with the full numerical calculations, but the 
other peaks are still red-shifted as compared with those 
found in the full numerical calculations. Finally, we in- 
clude one more level, i.e., |3) = 4 1 i ? , in the dynamics 
and hence obtain a four-level model. The results from 
such a four-level model are found to agree well with the 
full numerical results for the frequency region considered 
here. Detailed comparisons between these models are 
presented in Fig. [TJ 

The above-mentioned quantitative differences between 
various models can be explained by analyzing the transi- 
tion pathways in multi-level situations. Consider Eq. (ff| 
again. Now the term T' 3 ' for the three-level model in- 
cludes two integrals: 
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where the two integrals correspond to the following two 
sub-paths: |0) ->■ |1) ->■ |0) |1) (0-1-0-1) and |0) -> 
|1) |2> ->■ |1) (0-1-2-1). Clearly both sub-paths can 
contribute significantly to the total quantum transition 
amplitude from state |0) to state |1). Furthermore, the 
transition dipole moment //12 between states |1) and |2) 
is larger than the transition dipole moment /ioi between 
states |0) and |1) (see Table I). This indicates that the 
transition amplitude from |1) to |2) is more substantial 
than that from |0) to |1). It can be estimated that the 
contribution by the sub-path 0-1-2-1 of is about one 
order of magnitude larger than that by the sub-path 0- 
1-0-1. A detailed numerical comparison between them 
is shown in Fig. [5] Therefore the role of the sub-path 
0-1-0-1 in the previous two- level model should be re- 
placed by the sub-path of 0-1-2-1 in the present three- 
level model, thus explaining the quantitative difference 
between a two-level model and a three-level model. 

Similarly, for the interference between the pathways 
of and that accounts for the second peak at 
uj = 0.15 in Fig. [Ha), the contribution of the sub-path 
0-1-2-3-2-1 to the magnitude of is larger than the 
contributions of the sub-paths 0-1-2-1-2-1, 0-1-2-1-0-1, 
and 0-1-0-1-2-1. This can be interpreted by a four- level 
model. In addition, as for interference between multi- 
photon transitions of even higher orders, the results of 
the four-level model also agree well with the full numeri- 
cal calculations (Fig. 1(a)). This is because that the pop- 
ulations of higher states are so small under the present 
laser conditions that the additional sub-paths, including 
these higher states, are not important regardless of their 
larger transition dipole moments. 

Table I: Transition dipole moment of helium in atomic units. 
As discussed in the text, the big differences in the transition 
dipole moments are an important factor when explaining the 
quantitative differences between two-level, three-level, and 
four-level models. 





Ml2 


fJ-23 


0.4208 


2.499 


5.175 



Finally, we return to Fig. [IJb), where M. versus the 
laser's carrier frequency is shown for the Z X D state of he- 
lium. Similar to the case of 2 1 P, the peaks in Fig. QJb) 
can be associated with the interference between two 
neighboring transition paths at w = A_E/(2n), where 



I 

n = 1,2,..., and AE is the energy spacing between the 
ground state and the 3 1 !) state. Our general picture 
above can equally be applied to analyze this case and 
the details will not be repeated here. 



SUMMARY 

In summary, we have investigated numerically the fre- 
quency dependence of the carrier-envelope phase effect on 
bound-bound transitions of helium in an ultrashort laser 
pulse. It has been found that the CEP effect can occur 
regularly even at frequencies which are far off-resonance. 
To explain this numerical finding, we have examined a 
two-level model and developed a general and simple pic- 
ture, where the total transition amplitude can be decom- 
posed into different transition pathways. All these transi- 
tion pathways can be characterized by two indices n and 
j, with n being the number of photons involved and j 
indicating where the resonance frequency is located. For 
a long laser pulse, each pathway is associated with a nar- 
row peak at the resonance frequency LUj and the pathways 
at the same resonance frequency have the same depen- 
dence on the CEP. As a result, no interference can occur 
between different pathways and there is no CEP effect. 
In contrast, for an ultrashort pulse, the peaks for path- 
ways at different resonance frequencies can be broadened 
to cross with each other. Therefore, interference between 
different pathways can happen and lead to strong CEP 
effect. This general picture is valid for a wide range of 
laser intensities as long as the perturbation method is 
applicable, and can be generalized to multi-level models. 
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